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Abstract 

Nowadays, automotive engine compensation control is 
carried out electronically by utilizing many compensation 
maps in engine management systems; such that the engine 
can sustain its performance under the variations in engine 
operating conditions and environmental parameters. In 
traditional engine compensation map calibration, the 
parameters are normally set by a trial and error method 
because the exact mathematical engine model has not been 
derived. In this paper, a new framework, namely multi- 
input/output least-squares support vector committee 
machine, is proposed to construct the engine compensation 
control system (ECCS) models based on experimental data. 
As the number of adjustable parameters involved in the 
ECCS is very huge, the model accuracy and training time are 
usually degraded. Nonlinear regression is therefore 
employed to perform dimension reduction before modelling. 
The ECCS models are then embedded in an objective 
function for parameter optimization. Two widely-used 
evolutionary optimization algorithms, Genetic algorithm 
(GA) and particle swarm optimization (PSO), are applied to 
the objective function to determine the optimal calibration 
maps automatically. Experimental results show that the 
proposed modelling and optimization framework is effective 
and PSO is superior to the GA in compensation map 
calibration. 
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Introduction 

Nowadays, the automotive engine is controlled by the 
engine management system (EMS). The EMS 
calibration is achieved by parameterizing EMS 
software and looking for an optimal engine 
performance. The parameters stored in the EMS are 



organized in a map format. Calibrated maps under 
normal operating temperatures and conditions are 
called base maps. The base maps at least include base 
fuel map and base ignition map. Although the base 
map represents the fine-tuned performance of a 
vehicle engine, its performance also changes due to 
different engine operating and environmental 
conditions, such as barometric pressure (BP), engine 
temperature (ET), inlet air temperature (IAT) and 
battery voltage (BV). 

In order to maintain the engine performance under 
different environmental and engine operating 
conditions, every EMS has to include an engine 
compensation control system (ECCS) which usually 
organizes all control parameters in many 
compensation maps. The ECCS picks the trim values 
from the compensation maps based on the signals 
reported from the IAT, ET, knock and oxygen sensors. 
The trim values are then added to or deducted from 
the fuel injection duration and ignition timing of the 
base fuel and ignition maps to become the final engine 
control output. This process provides an adaptive 
engine control action. In addition, the ECCS should 
essentially allow the engine to achieve rapid warm-up, 
good acceleration and deceleration responses, fair 
knock intensity, short cranking time (CT), sustainable 
air-ratio (lambda value) and torque under various 
weather and operating conditions. 

In engine compensation control (ECC), a 
compensation map is represented as a row of trim 
values that can normally be visualized in form of 
curves. There are many compensation maps in an 
ECCS, including (1) the fuel compensation maps for 
IAT, ET, BP, BV and cranking; (2) the ignition 
compensation maps for IAT, ET. An example of fuel 
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FIG. 1 AN EXAMPLE OF FUEL COMPENSATION MAP FOR ET WHERE THE ET IS DISCRETELY DIVIDED 



compensation map and curve for ET is shown in Fig. 1 
(MoTeC, 2002); where the engine temperature is 
discretely divided into o C,l0 o C,...,ll0 o C, and the values 
in the table represent the percentages of fuel trims, 
which are tunable parameters. In practice, all 
compensation curves or trim values are engine 
dependent. Every car model has to undergo a 
complicated calibration process of compensation maps 
during development. 

The ultimate goal of the engine compensation map 
calibration (ECMC) is to tune up the engine optimally 
under different environmental and engine conditions 
but the goal is extremely difficult to achieve. Current 
practice of ECMC relies on the experience of the 
automotive engineer who handles a huge number of 
combinations of engine control parameters. The 
relationship between the input and output variables of 
a modern car engine is a complicated multi-variable 
nonlinear function (Robert Bosch, 2004; Li, 2007) 
because the automotive engine is an integrated system 
of thermo-fluid, electromechanical and computer 
control systems. Consequently, ECMC is usually 
performed empirically by the automotive engineer 
through many experiments (Robert Bosch, 2004; Li, 
2007). Therefore, vehicle manufacturers normally 
spend up to 12 months to tune up an EMS optimally 
for a new car model (Bell, 2006). 

To date, a few studies on computer-aided calibration 
for some simple EMS parameters are available. 
However, research on computer-aided calibration for 
engine compensation maps is hardly found in the 
open literature. The most common calibration 
approach for simple EMS parameters is simulation- 
based calibration (Rask & Sellnau, 2004; Tan et al, 
2006)), which divides the calibration process into two 
parts: modeling and optimization. In terms of 
modeling, a mathematical engine model can be 
constructed either from empirical models (Saerens et 
al, 2009; Zhang et al, 2010; Errico et al, 2011; Qi et al, 
2011) or commercial engine simulation software such 



as Engine Analyzer Pro, GT Power and WAVE'S 1-D 
cycle simulation, etc (Wu et al, 2007). The empirical 
models and simulation software are derived from 
some simple physical laws. These mathematical 
models can simulate engine performance 
approximately, and thus can be employed as objective 
functions. The optimal EMS parameters can then be 
determined based on the objective functions by using 
computer-aided optimization methods such as 
gradient search methods (Wu et al, 2007) and genetic 
algorithm (GA) (Errico et al, 2011). Hence, the amount 
of time and resources required for engine 
development can be significantly reduced. However, 
the predictive accuracy of the mathematical engine 
models is poor due to the following facts: 

1. To predict engine performance, the simulation 
software and mathematical models require the 
engineer to provide many engine-specific 
parameters such as combustion efficiency. In 
practice, it is too demanding to obtain or estimate 
such parameters. 

2. Excessive numbers of assumptions have been 
made in these models. This renders such engine 
models overly simple in comparison with real 
engine systems. 

Several recent studies have described the use of neural 
networks (NNs) to process collected experimental data 
and generate engine performance models (Celik & 
Arcaklioglu, 2005; Garcia-Nieto et al, 2008; Dickinson 
and Shenton, 2009; Togun and Baysec, 2010). In this 
way, many direct search techniques can then be 
applied to search for the optimal engine parameters 
automatically. However, there are two main 
drawbacks to the NN modeling approach: 

1. The architecture, including the number of hidden 
neurons, must be determined a priori or modified 
during training by heuristics, which results in a 
sub-optimal network structure. 

2. During the training process, NNs can easily 
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become stuck in local minima. Various methods 
of avoiding local minima (e.g. early stopping and 
weight decay) have been employed. However, 
these methods greatly affect the generalization of 
the estimated function (Haykin, 1999). 

Recently, an emerging technique of least-squares 
support vector machines (LS-SVM) (Suykens et al, 
2002) has been applied to model the engine 
performance (Vong et al, 2006). LS-SVM combines the 
advantages of neural networks (particularly their 
ability to handle a large amount of highly nonlinear 
data) and nonlinear regression (with its high capacity 
for generalization). It can achieve better accuracy than 
traditional NNs for automotive engine performance 
modelling. Subsequently, the authors proposed the 
application of multi-input/output least-squares 
support vector machines (MIO LS-SVM) and a GA to 
determine the optimal set-points for a proportional- 
integral-derivative idle-air valve controller and some 
engine control variables for automotive engine idle- 
speed control (Wong, Tarn et al, 2008), which is a very 
difficult multi-objective optimization problem. The 
integrated approach can overcome all of the 
limitations of the existing EMS calibration methods 
mentioned above. However it is not suitable for 
modelling the ECCS because the input variables of the 
ECCS model are usually a mixture of tunable and 
untunable parameters. The untunable parameters are 
uncontrollable environmental and operating factors 
such as ET, IAT, BV, and BP, which cannot be directly 
regulated by the optimization algorithm. Nevertheless, 
these untunable parameters must be taken into 
account for obtaining optimal engine performance at 
different environmental and engine operating 
conditions. As the MIO LS-SVM+GA methods (Wong, 
Tarn et al, 2008) can deal with tunable parameters only, 
this paper proposes a new framework utilizing MIO 
LS-SVM and committee machine (Krogh & Vedelsby, 
1995; Haykin 1999; Chen et al, 2009), namely, multi- 
input/output least-squares support vector committee 
machine (MIO LS-SVCM) to model the complex ECCS 
as an objective function, and then proper optimizer 
can perform optimization based on this kind of 
objective function. 

Proposed Modelling and Optimization 
Framework 

The MIO LS-SVCM framework is inspired from the 
concept of committee machines that a group of experts 
(models) can produce a better result relative to a single 
expert; especially when the training data of an 



application can be separated independently. However, 
there is no universal formulation for committee 
machines, the application of committee machines is 
case dependent (Haykin, 1999). For this application, 
the detail formulation is specially designed. The 
proposed MIO LS-SVCM firstly separates the whole 
ECCS model into some simpler sub-models according 
to the feasible combinations of untunable parameters, 
such as IAT, ET and BV. In this treatment, the sub- 
models do not contain any untunable parameters that 
enable the compensation map optimization. In order 
to restrict the number of combinations of the 
untunable parameters, the values of the untunable 
parameters should be bounded and discretely divided. 
This idea matches with the general compensation map 
format of the EMS (Fig. 1). As there are still many 
adjustable parameters for the divided ECC sub-models, 
the model accuracy and training time are usually 
degraded. However, the relationships of adjustable 
parameters in some compensation maps could be 
expressed as compensation curves, so nonlinear 
regression can be applied to represent the data points 
of curves by several coefficients so as to realize 
dimension reduction before modelling. Afterwards, 
the multi-input/output LS-SVM can be applied to 
construct the ECC sub-models individually. Each ECC 
sub-model represents one engine output performance 
at a specific combination of engine operating and 
environmental condition. Finally, the ECC sub-models 
are weighted and systemically aggregated to form an 
objective function (i.e. committee machine) for ECMC. 

A schematic illustration of the framework and overall 
methodology of LS-SVCM is shown in Fig. 2. The 
upper branch in Fig. 2 shows the steps required to 
construct the MIO LS-SVCM model. Experiments are 
still required but only to provide sufficient data for LS- 
SVM training and validation. The design of 
experiments (DOE) methodology is additionally used 
to streamline the process of creating representative 
sampling data points to train the models. 
Normalization is applied to preprocess the tunable 
parameters in the experimental datasets. Once the 
ECC sub-models obtained are validated in terms of 
accuracy, it becomes possible to use a proper 
computer-aided optimization technique to search for 
the best engine calibration maps automatically. The 
entire model derived by MIO LS-SVCM is obviously 
discontinuous such that it is non-differentiable, and 
gradient information of this model cannot be obtained. 
In practice, evolutionary optimization methods are 
usually employed to determine the optimal solution 
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FIG. 2 FRAMEWORK FOR MODELING AND OPTIMIZATION OF ECC PARAMETERS 



for this kind of discontinuous problem because no 
gradient information of the model is required. 
Evolutionary optimization in general is methods that 
simulate natural evolution for the task of global 
optimization. However, there are many evolutionary 
optimization algorithms available in the open 
literature, it is impossible to exam all of them. In order 
to check if the proposed framework can fit various 
evolutionary optimization algorithms, two widely- 
used evolutionary optimization techniques, genetic 
algorithm (GA) and particle swarm optimization 
(PSO), were selected and tested for this large-scale 
constrained multivariable optimization problem. GA 
has already been a well-known technique for solving 
many automotive engineering optimization problems 
(Chiou & Liu, 2009; Gauchia et al., 2010). PSO 
(Kennedy & Eberhart, 1995) is also an evolutionary 
computation technique similar to the GA. However, it 
has no evolutionary operators of crossover and 
mutation as with the GA. From the user point of view, 
PSO can reduce the time and burden of trying 
different evolutionary operators to find the optimal 
solution. 



In the lower branch of Fig. 2, the optimal ECC set- 
points by both optimizers are denormalized into true 
values and remapped into compensation maps 
according to the format of the EMS. Finally, the true 
set-points in the maps are loaded to the EMS to carry 
out evaluation tests and comparison. The feasibility 
and efficiency of the two optimization approaches can 
therefore be examined. The purpose of this paper is to 
demonstrate the effectiveness of the proposed MIO 
LS-SVCM + GA/PSO methods on ECMC problems. 
The detail operation of the modelling and 
optimization framework and the specific techniques 
are described in the following sections. 

Engine Compensation Control Model 
Identification and Case Study 

LS-SVM Formulation for Multi-class and 
Multivariable Function Construction 

The typical LS-SVM modelling algorithm is a multi- 
input and single-output modelling method. To deal 
with the current ECMC problem, the MIO LS-SVCM 
model is constructed based on many typical LS-SVM 
models. The concept is presented below. 
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FIG. 3 MIO LS-SVCM MODELLING FRAMEWORK 



Consider a training dataset, D = {(x k ,y k j )}^, ! : =l , with 

N data points where Xk e R n represents the k"> ECC 
parameter set-up (i.e. set-points), and yk,i e R m is the k th 
corresponding engine output performance, k = 1 to N 
and i = 1 to z. The index represents a reasonable 
combintation of untunable environmental parameters 
or engine operating conditions. In this study, this set 
of untunable parameters contains the engine operating 
conditions of ET, IAT, and BV only. For example, 
when i = 1, the combination of untunable parameters is 
ET = 20°C, IAT = 20°C, and BV = 12V. Here, y*,; is an tri- 



dimensional output performance vector, yk,,= 
[yk,i,i,...,yk, m ,i]. In this applicaton, yu,; could be the engine 
knock intensity at the i th engine operating condition 
and yk,m,i could be the cranking time at the i th engine 
operating condition. For an automotive engine, each 
engine output performance in yi,i is usually an 
individual variable and can be measured separately, 
so the training dataset D can be arranged 
independently as follows: 

D = |(Du,... D M ,...Dm,i),.. (Dy,... Da,/,... Dm,/),.. (Du,... 
D>,,z,...D ra ,z)}, (1) 
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where D h . = {{x k ,y kXi )} N k=l , z i=l ; h e [l,m]. In other words, 
the entire training dataset D is separated into z classes 
x m multi-input/single-output sub-training datasets 
Dh,i (Fig. 3). For each multi-input/single-output dataset 
Um, the LS-SVM dual formulation for nonlinear 
function estimation as given in Equation (2) (Suykens 
et al., 2002) can be applied to. 



Solve in a h .,b h y. 
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Where In is an N-dimensional identity matrix, yi,,i = 
[y\,Ki, yN,h,i] T , It. is an N-dimensional vector = [1,...,1] T , 
and ach,i =[au,i,-, aN,h,i] T . The resulting LS-SVM model 
for function estimation becomes: 



M Ki (x) = ^a kXl <f>(x k ) T <f>(x) + b hl 

k = \ 

N 

= Y, a kxMx k ,x) + b hi 
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where aw, hjeR are the solutions of Equation (2), x is 
the new input ECC set-points for engine performance 
prediction. After inferring m x z pairs of 
hyperparameters cn>,<) by Bayesian inference (Vong 
et al., 2006), m x z individual training datasets are used 
to calculate m x z individual sets of support vectors 
afc/i,;and threshold values bhj. Finally, m * z individual 
sets of Mh,i(x) are constructed based on Equation (3). 
Each Mh,i(x) sub-model is viewed as a committe 
member in the committee machine. The whole MIO 
LS-SVCM modelling algorithm is shown in Fig. 3, 
where a set of typical LS-SVM sub-models are 
generated to predict the engine output performances 
at the i th reasonable combination of engine operating 
condition for different combinations of ECC variables 
in x. 

Experimental Set-up for Case Study 

In this case study, a Honda Integra DC5 Type-R with a 
K20A DOHC i-VTEC engine was used as a test car to 
demonstrate the feasibility of the proposed 
methodology through a programmable MoTeC M800 
EMS (Fig. 4). As an illustration, six compensation 
maps were considered in the experiment, including 
the fuel compensation maps for IAT, ET, CT and BV, 



as well as ignition compensation maps for IAT and ET. 
The fuel compensation map for BP was not considered 
because the base map of the test engine uses speed- 
density system, which estimates engine load by 
measuring manifold-absolute pressure. The speed- 
density system can therefore automatically correct for 
BP changes (Hartman, 1993). Moreover, four engine 
output performance data were considered including 
knock intensity (knk), cranking time, differential 
lambda (AA) and differential engine torque (AT). The 
lambda value and knock intensity were measured by a 
wide-band lambda sensor and a knock monitor 
respectively. The cranking time was recorded by the 
built-in data logging function in the MoTeC EMS. The 
differential engine torque was measured by running 
the test car on a chassis dynamometer (dyno) at an idle 
speed and full-load condition at random engine 
speeds, so that various conditions of ET and IAT were 
achieved for data sampling. All of the experimental 
data were collected under a controlled ambient 
environment using powerful air-conditioners and fans. 

Selection of Ranges 

In this case study, the fuel and ignition compensation 
maps are the functions of the untunable 
environmental parameters IAT, ET, and BV. Owing to 
the limitations in the automotive laboratory 
equipment at the University of Macau, test car 
characteristics, and weather in Macau, the sampling 
temperature and BV ranges were selected as: 15 °C to 
55°C for IAT, 15°C to 105°C for ET and 11.5V to 13.5V 
for BV. Based on the format of the programmable EMS, 
the intervals for IAT and ET were set at 10 °C, i.e. IAT= 
(20, 30,... ,50), ET = (20, 30,... ,100) and for BV, the 
interval was set at 1, i.e. BV = (12, 13). 

According to the concept of MIO LS-SVCM, the 
complete ECCS model is therefore separated into sub- 
models according to different temperature and BV 
intervals. The indexes (e, a, b) were used directly as 
indexes for ET, IAT and BV respectively. The indexes 
of e, a, and b covered the range for e ± 5°C, a ± 5°C and 
b ± 0.5V respectively. For example, the data from 15 °C 
to 25 °C were viewed as 20 °C, then 20 was used as one 
of the model indexes. 

Originally, the number of combinations of such 
indexes is 72. However, many combinations of IATs 
and ETs are unreasonable. For example, it is 
impossible to have 20 °C IAT and 100 °C ET 
simultaneously. By considering the real situations, the 
number of feasible combinations was reduced to 30 (i.e. 
15 feasible combinations of IATs and ETs x 2 BV 
combinations) as shown in Table 1. 
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FIG. 4 EXPERIMENTAL SET-UP FOR DATA SAMPLING AND PROGRAMMABLE EMS 
TABLE 1 FEASIBLE COMBINATIONS OF IATs AND ETs WITH BVs = 12V AND 13V 
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(1: feasible; 0: unreasonable) 
TABLE 2 FEASIBLE INPUT VARIABLES FOR THE USER TO FILL IN THE COMPENSATION MAPS 
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Dimension Reduction of Input Variables by Nonlinear 
Regression 

According to the selected maps, ranges, intervals, 
indexes and the format of the EMS, there were six 
compensation maps with 39 input variables (Table 2). 
However, 39 variables are still too much. To avoid 
handling more of such variables, each fuel and 
ignition map for IAT and ET were converted into 
compensation curves with 4 coefficients, which 
uu, Cu f du with u =1,2,3,4 , using nonlinear regression. 
For each curve, the x-axis represents the value of IAT 
and ET, and the y-axis represents the value of 
compensation factors of fuel or ignition. The 



compensation factors for fuel and ignition are 
percentages of fuel trim and degrees of ignition angle 
trim respectively. One more digit, p u e{0, 1}; u e [1,4], 
was used as identity digit to represent the categories 
of curves for IAT and ET. By referring to (Hartman, 
1993; MoTeC, 2002; Wang et al„ 2006, Wong, Vong et 
ah, 2008), the compensation curves for cranking and 
BV can also be expressed in linear and hyperbolic 
functions respectively. Therefore, the input x was 
reduced from 39 to 23 as shown in Equation (4). The 
parameters in the vector x are tunable/adjustable 
parameters. 

x = < ai, bi, a, di, pi, ai, bi, a, di, pi, en, fo, a, di, pi, cu, bi, 
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a, di, pi, as, a6, be > (4) 
where 

ai, bi, a, dr. Coefficients of fuel curve for IAT 
compensation 

pr. Identity digit of fuel curve for IAT compensation 

02, bi, C2, dr. Coefficients of fuel curve for ET 
compensation 

pr. Identity digit of fuel curve for ET compensation 

as, bi, a, dr. Coefficients of ignition curve for IAT 
compensation 

pr. Identity digit of ignition curve for IAT 
compensation 

ai, bi, a, dr. Coefficients of ignition curve for ET 
compensation 

pr. Identity digit of ignition curve for ET compensation 
as: Coefficients of fuel line for cranking compensation 
as, be: Coefficients of fuel curve for BV compensation 
Overall, the shape of the compensation curves could 



generally be represented as a third order polynomial, 
Equation (5), or sinusoidal function, Equation (6), or 
hyperbolic function, Equation (7) according to many 
text books and user manuals (Hartman, 1993; MoTeC, 
2002; Wang et ah, 2006). Linear regression can also be 
obtained using the polynomial regression of Equation 
(5) by setting the first two coefficients to zero. 



Y = aX+b u X 2 +c u X+d u 

Y = a u +b u (cos(c u X) + d u ) 



Y = *- + b K 



X 



(5) 
(6) 

(7) 



Where X and Y are the map indexes (e.g. at 30 °C ET, X 
= 30) and compensation factor (e.g. for -4.5% of fuel 
trim, Y = -4.5) in the compensation maps respectively. 
In the input vector x, Equations (5) and (6) are 
respectively represented by setting p u to be and 1 for 
each u. The use of the digital number as an extra 
tunable parameter allows the computer optimizers to 
evaluate and determine the best shape of 
compensation curve in the course of optimization. 



TABLE 3 DIVISION OF DATASET D INTO 30 SUB-DATASETS Dw ACCORDING TO VARIOUS ETs, IATs and BVs 
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-0.2587 





0.2345 


15.9 


5.3 


0.078 


3.1 


30 


1 


D H.100S0M 


5.1428 


-0.2587 





0.2345 


15.7 


2.7 


0.088 


-10.7 


1 


2 


D )i,20,20,12 


7.7142 


-0.4607 


1 


0.2675 


15.8 


4.9 


0.078 


3.1 


2 


2 


D M0,20,12 


7.7142 


-0.4607 


1 


0.2675 


15.8 


3.4 


0.225 


-3.2 


30 


2 


D 1,100,50,13 


7.7142 


-0.4607 


1 


0.2675 


15.6 


3.4 


0.027 


1 


1 


300 


D h, 20,20,12 


9.1428 


-0.4964 


1 


0.2126 


15.8 


3.3 


0.225 


3.2 


2 


300 


D fc,30,20,12 


9.1428 


-0.4964 


1 


0.2126 


16.4 


3.2 


0.153 


-3.5 


30 


300 


D (1,100,50,13 


9.1428 


-0.4964 


1 


0.2126 


16.1 


1.6 


0.09 


-8.4 
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Data Sampling and Division for Case Study 

In view of the structure of MIO LS-SVCM, the data 
was sampled and divided according to different 
feasible combinations of IATs, ETs and BVs. Moreover, 
the training dataset was expressed as: D*,FDi,e^=(«, 
yKi \=\xk„ yk,e,a,b}= {x k„ y k,h,e,a,b\, where k = 1 to 300 in this 
case study. The index i represents the i th combination 
of (e, a, b). yk.h.e.a.b is the h lh engine performance data in 
the k th sample data point at e, a and b. Specifically, yn is 
defined as follows: 

yk,i= {knh, CTi, AAi, AT,) = \knke,a,b, CTe,a,b, AAe.a.b, ATe,a,b} 

(8) 

and 

A^U-A'I (9) 

AT.= Ti-T (10) 

For the above datasets, knh is the engine knock 
intensity, CTi is the engine cranking time, AA; is the 
error between the lambda value produced by the base 
maps and the lambda value after compensation, AT/ is 
the error between the output torque produced by the 
base maps and the torque after compensation. A 'and 
T are respectively the lambda value and the output 
torque at normal operating temperatures e = 90 °C, a = 
40 °C and b = 13V without compensation. According to 
Table 1, the feasible combination of the training data 
and LS-SVM models are shown below: 

i th combination = (e,a,b) = 

[20,20,12 30,20,12 40,20,12 50,20,12 60,20,12 70,20,12 
50,30,12 60,30,12 70,30,12 80,30,12 90,30,12 100,30,12 

90.40.12 100,40,12 100,50,12 20,20,13 30,20,13 40,20,13 

50.20.13 60,20,13 70,20,13 50,30,13 60,30,13 70,30,13 

80,30,13 90,30,13 100,30,13 90,40,13 100,40,13 100,50,13] 

(11) 

In the above representation, when z'=l, (e,a,b) = 

(20,20,12); z=2, (e,a,b) = (30,20,12); ; z=30, (e,a,b) = 

(100,50,13). According to above combination, D could 
be grouped into 30 training model datasets, namely 
Dfc,2o,2o,j2, D h,3o,2o,u, T) h,wo,5o,i3 . The training dataset D 
is expressed in Table 3. 

Latin hypercube sampling (Lunani et al, 1995) is a 
technique of DOE which was employed to select a 
representative set of operating points to generate 
training samples. A total of 300 sets of representative 
combinations of input variables were selected and 
loaded to the EMS to produce 36000 sets (i.e. 300 data 
x 30 training model datasets * 4 output performance 



values) of engine output performance data. After 
collecting a sample dataset D, every data subset Dwc 
D was randomly divided into two sets: TRAIN/,,, for 
training and TEST/,,,- for testing, such that Dw = 
TRAINwUTESTm, where TRAIN/.,; contained 80% of D«, 
TEST/,,, held the remaining 20% for h =[1,4], and the 
engine performance for each h in this case study is 
shown in Table 3. Subsequently every TRAIN/,,, was 
sent to the LS-SVM program for training which was 
implemented using MATLAB on a PC. 

Data Normalization 

To obtain a more accurate modelling result, the input 
data of the training dataset is conventionally 
normalized before training (Pyle, 1999). This 
procedure prevents any input parameter from 
dominating the output value. For every input 
parameter, it was normalized to within the range [0, 1], 
i.e., a unit variance, using Equation (12). 



v — V 
max min 

The limits of Vxmn and i/max for each variable were 
obtained from handbooks or expert knowledge. After 
obtaining the optimal setting, each set-point must be 
de-normalized to obtain the actual value v. The 
processes of normalization (v — > v") and de- 
normalization (c* — > v) were performed by using 
Equation (12). The process flow of the normalization 
and de-normalization is shown in Fig. 2. 

Modelling Results of MIO LS-SVCM 

After data division and pre-processing, each subset 
Dh,i was passed to the LS-SVM program in Equation (2) 
in order to construct the function Mh,i(x) in Equation 
(3). According to the division of the training data, 
there were totally 120 output models trained for every 
h and i (i.e. 4 performance outputs x 30 classes of 
models). In other words, the LS-SVM program was 
run for 120 times. In every run, a different subset 
TRAIN/,,, was used as the training dataset to construct 
the corresponding engine output performance 
function. The one hundred and twenty trained engine 
performance functions or models Mh,i(x) are grouped 
as follows: 

[Mh,i(x)] = [Mh,e,a,b(x)] =[Mh,20,20,12(x) Mh,30,20,u(x) Mh,100,50,13(x)] 

(13) 

Where h is the index of output performance function 
of the knock intensity, cranking time, air-ratio 
difference and engine torque difference, and they are 
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given in Equations (14) to (17). 
For h=l, Mi,i(x) = knk(x) (14) 
For h=2, M2,i(x) = CT(x) (15) 
For h=3, Mxi{x) =AAj(x) = I Ai(x) - X' I (16) 
For h=4, M4,i(x) =ATi(x)= Ti(x)-T (17) 

where 

knk(x)= knke, a ,b(x): Knock intensity function for engine 
temperature e, inlet air temperature a and battery 
voltage b. 

CTi(x) = CTe,a,b(x): Cranking time function for e, a and b. 

Ai(x) = Ae A ,b(x): Lambda function for e, a and b. 

Ti(x) = Te,a,b(x): Output torque function for e, a and b. 

In the aforementioned functions, e,a,b are used as a 
function index again. For example: when Ae,a,b(x) = A so, 
30, u(x), it represents a lambda function for the engine 



TABLE 4 ACCURACIES OF DIFFERENT ENGINE OUTPUT PERFORMANCE FUNCTIONS 



i 


Engine output 
function index 

e,a,b 


Mean accuracy of 
knh withTESTi,, (%) 


Mean accuracy of 
CTi withTEST2,,(%) 


Mean accuracy of AA< 
withTEST3,,(%) 


Mean accuracy of 
ATi with TEST*,; 
(%) 


1 


20,20,12 


98.97 


90.85 


92.24 


96.77 


2 


30,20,12 


99.06 


83.36 


94.74 


96.44 


3 


40,20,12 


98.00 


92.44 


93.36 


90.19 


4 


50,20,12 


99.04 


88.85 


98.44 


94.17 


5 


60,20,12 


98.82 


92.87 


97.74 


92.45 


6 


70,20,12 


97.49 


88.54 


93.72 


93.72 


7 


50,30,12 


98.08 


91.31 


94.25 


95.30 


8 


60,30,12 


98.56 


83.03 


90.24 


96.06 


9 


70,30,12 


99.24 


85.73 


97.20 


94.87 


10 


80,30,12 


98.56 


91.09 


91.53 


92.81 


11 


90,30,12 


97.23 


92.73 


94.87 


91.24 


12 


100,30,12 


98.67 


91.35 


92.81 


93.42 


13 


90,40,12 


98.34 


86.93 


91.24 


96.67 


14 


100,40,12 


97.64 


84.38 


93.42 


97.57 


15 


100,50,12 


99.08 


86.47 


96.67 


95.89 


16 


20,20,13 


98.37 


91.25 


97.57 


96.63 


17 


30,20,13 


98.96 


92.70 


95.89 


93.51 


18 


40,20,13 


98.43 


86.24 


96.63 


95.44 


19 


50,20,13 


97.24 


83.72 


93.51 


98.36 


20 


60,20,13 


98.82 


84.90 


95.44 


96.34 


21 


70,20,13 


97.69 


87.45 


98.36 


97.74 


22 


50,30,13 


98.38 


93.61 


96.34 


93.72 


23 


60,30,13 


98.56 


90.46 


95.25 


94.25 


24 


70,30,13 


99.24 


93.35 


93.14 


90.24 


25 


80,30,13 


98.56 


85.32 


99.34 


97.20 


26 


90,30,13 


98.23 


84.89 


94.87 


91.53 


27 


100,30,13 


98.67 


91.63 


92.81 


94.87 


28 


90,40,13 


98.34 


83.65 


91.24 


92.81 


29 


100,40,13 


97.64 


93.12 


93.42 


91.24 


30 


100,50,13 


98.48 


90.03 


96.67 


93.42 


Average 




98.41 


88.74 


94.77 


94.50 



Overall average: 94.10% 



temperature e ~ 80 °C, inlet air temperature a ~ 30 °C 
and battery voltage b ~ 12V. The accuracy of each 
function of Mw (x) was verified using an error function 
proposed in Equation (18). 




(18) 



where xt is the ECC parameter set-up of the t th data 
point in TEST*,;, yt.hj is the actual engine output value 
of the f th data point {xt, yt.hj] in TEST*,,, and Nt is the 
number of data points in the test dataset. The error Eh,i 
is a root-mean-square of the difference between the 
true value yt,h,i and its corresponding value estimated 
from Mh,i(xt). The difference is also divided by yw, so 
that the result is normalized to within the range [0,1]- 
Hence, the accuracy rate for each output function 
Mhj(x) is calculated using Equation (19). 



Accuracyhj = (1 - Em) x100% 



(19) 
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According to the indexes i and h, each class of engine 
performance functions contains 7200 (i.e. 4 output 
performance values x 300 data x 30 classes of models x 
20% ) test data in TEST*,/. Based on the index h, each 
engine performance function uses 1800 (i.e. 7200 data 
/4 output performance values) test data for evaluation. 
By comparing TESTw with the predicated output 
performance data from the corresponding engine 
performance function, the accuracy for each function 
was obtained as shown in Tables 4. The results 
indicate that the overall function accuracy is very good, 
at 94.1%. However, it is believed that the accuracy of 
the function can be improved by increasing the 
amount of training data. 

Optimization of Compensation Maps 

After obtaining the ECCS models by MIO LS-SVM, the 
models were aggregated to form an objective function. 
It was then possible to use the GA and PSO to search 
for the best ECC parameter set-up automatically. In 
this application, the ECC set-up invloves 23 tunable 
parameters with multiple optimization objectives and 
constraints, so it is a challenging large-scale 
optimization problem. To minimize the chance of 
obtaining a local optimum, several different and 
arbitrarily generated initial populations are used. For 
this purpose, ten independent optimization runs were 
performed for each optimizaiton method. The results 
of both optimization algorithms are also compared at 
the end of this section. 

Objective Function for ECMC 

An objective function was systemically designed to 
evaluate the engine performance under different 
compensation control set-ups. In this case study, the 
main compensation objective is to minimize the 
following four basic requirements (Jurgen, 1995) and 
subject to some constraints: 

(1) Fair knock intensity 

Higher knock intensity means higher engine 
torque output. However, high knock intensity 
decreases the service life of an engine. Empirically, 
the knock intensity should be adjusted below 60% 
for common vehicles and the value can be 
adjusted by the user. 

(2) Short cranking time 

Under a fine-tuned base map, vehicles are 
expected to start as rapidly as possible in any 
weather conditions. 

(3) Sustainable air ratio (Ideal value: AA=0) 



The engine compensation control must be capable 
of maintaining the output lambda value close to 
the lambda value specified by the base maps 
under different engine operating conditions. 

(4) Sustainable engine torque (Ideal value: AT^O) 

This is a balance between the engine torque and 
knock intensity. When the ignition advance drops 
down, or becomes a zero value, the knock 
intensity decreases. However, the engine torque is 
also decreased significantly. Therefore, a function 
is required to avoid excessive ignition retard in the 
ECMC. 

According to the above requirements, the objective 
function for this case study is to minimize the sub- 
functions from the following equations: 



f a (x) = AA = 



30 



f b (x)= CT 



f c (x) = AT-. 



£ct; (x) 



30 



to- 7 ;') 



30 



(20) 



(21) 



(22) 



Hence, the overall objective function /<*; is defined as: 

fa = Max\fitness] =Max[-W r itan- 1 (f,(x))-W 2 tan- 1 (/i,(a:)) 

+W 3 tan<(fc(x))] (23) 

Wi,W2,W3 are the user-defined weights of fi(x), fi(x) 
and fc(x) respectively, where fa(x) is the average air- 
ratio difference, fi(x) is the average cranking time and 
fc(x) is the average engine torque difference, subject to 
the following constraint equations: 



Vknh(x)<60%,i=l,2, ,30 

V« e{0,l} 



(24) 
(25) 



Equation (25) ensures that each selected compensation 
curve by the optimizer is confined to the shapes of the 
curves described by Equations (5) to (6). In this case 
study, Wi, Wi and W3 were all set at 1. Each 
performance index was also transformed to a scale of 
(0, n/2) by tan-!(-) in Equation (23). This ensures that 
each index makes the same contribution to the 
objective function. Since the objective function /<%• 
considers many ECMC requirements and constraints, 
it is treated as the overall compensation performance 
of the EMS. 
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TABLE 5 GA OPERATORS AND PARAMETERS 



Number of generation 1000 
Population size 50 

Selection method Standard proportional selection 

Crossover method Simple crossover with probability = 80% 

Mutation method Hybrid static Gaussian and uniform mutation with probability 

= 40% and standard deviation = 0.2 



TABLE 6 GA OPTIMIZATION RESULTS OF DIFFERENT OUTPUT FUNCTIONS 



Engine output 


Knki(x)(%) 


CTi(x)(s) 


AAi(x) 


ATi(x)(Nm) 


function Mh,i(x) 


h=l 




h=3 


h=4 


Mh,20,20,u(x) 


1 £ no 
io.uy 


n ao 

U.4Z 


n ai 

U.41 


-U.1D 


Mh,30,20,u(x) 


1 £ AQ 
10.4y 


i ift 

J.JO 


n ai 

U.41 


1 nn 

l.UU 


Mh,40,20,u(x) 


lO.OO 


1 Q7 
1.(5/ 


n 10. 
u.oy 


O.Uo 


Mh,50,20,u(x) 


1 £ 01 
lO.Zl 


1 A£ 
1.40 


n o^ 

U.ZO 


o nn 
z.UU 


M.h,60,20,u(x) 


10. ZZ 


o no 
z.Uo 


n 17 
U.O/ 


1 QQ 

i.yy 


M.h,70,20,u(x) 


1 71 
10./ 1 


1 Q1 


n iq 
u.oy 


n ir 
U.oo 


Mh,50,30,12(x) 


1 £ 10 

10. Oz 


£0 
Z.Dz 


n 17 
u.o/ 


n rs3 
U.Oo 


M.h,60,30,12(x) 


1 £ 10 
ID. OZ 


1 A1 
1.41 


n is 

U.Oo 


n n^ 
-u.uo 


M h,70,30,l2(x) 


1 £ 10 
ID. Oz 


111 
0.10 


n 7Q 
u./y 


n 07 
u.z/ 


Mh,80,30,u(x) 


1 c no. 
io.uy 


O 

Z.00 


u.oo 


1 77 
1.// 


Mh,90,30,u(x) 


1D.D/ 


1 77 


111 
i.ii 


1 

i.oy 


Mh,100,30,12(x) 


1 £ £D 
1D.DU 


117 
1.1/ 


n so 

U.oz 


n n^ 
-u.uo 


Mh,90,40,12(x) 


1 A 70 
10. /Z 


o no 
z.uz 


U.jO 


n ^ 
-u.oo 


Mh,W0,40,12(x) 


10.1/ 


O 7Q 

z./o 


1 OQ 
l.ZO 


n 77 
U./ / 


Mh,100,50,12(x) 


ID. 04 


1 1/1 
1.04 


n /i7 
U.4/ 


n ic 
U.OO 


Mh,20,20,13(x) 


1 £ 1 Q 

id. iy 


U.40 


n /ii 
U.40 


n 01 
-U.zl 


A/I it 3(1 7(1 1 2 ( Y ) 


16.27 


3.38 


0.43 


0.80 


M.k,40,20,is(x) 


16.46 


1.87 


0.41 


6.28 


Mh,50,20,13(x) 


16.31 


1.46 


0.28 


2.20 


Mh,60,20,13(x) 


16.32 


1.08 


0.39 


1.79 


Mh,70,20,13(x) 


16.81 


0.83 


0.41 


0.55 


Mh,50,30,13(x) 


16.42 


1.62 


0.39 


0.78 


Mh,60,30,13(x) 


16.42 


1.41 


0.41 


0.14 


Mh,70,30,13(x) 


16.42 


1.13 


0.82 


-0.47 


Mh,80,30,13(x) 


16.03 


2.11 


0.64 


1.56 


Mh,90,30,13(x) 


16.45 


1.88 


0.47 


1.79 


Mh r 100,30,13(x) 


16.17 


3.38 


0.43 


-1.00 


Mh,90,4O,13(x) 


16.82 


1.02 


0.42 


-0.12 


Mh,100,40,13(x) 


17.03 


0.51 


0.45 


-0.12 


Mh,100,50,13(x) 


16.64 


0.34 


0.49 


0.55 


Overall average 


16.42 


1.75 


0.52 


1.02 


Fitness value 






-39.5 
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GA Setting and Optimization Results 

The GA optimization framework was implemented 
using MATLAB. By testing various crossover and 
mutation methods for this application, the GA 
operators and parameters as shown in Table 5 were 
selected to ensure maximum efficiency and accuracy. 
Ten sets of optimization results were obtained by the 
GA, which were different from each other. The 
standard deviation of the fitness of the ten 
independent runs is about 0.06. Table 6 shows that the 
best fitness value out of the ten independent runs is - 
39.5. The optimized set-points recommended by the 
GA are then de-normalized into true values and 
remapped into six compensation maps according to 
the format of the EMS and Equations (5) to (7). The 
optimal maps are shown in Table 7. 

PSO Setting and Results 

The PSO framework was again implemented using 
MATLAB. In PSO, the population is called a swarm, 
and the individuals (i.e. different combinations of ECC 
set-points x) are called particles. The selection of the 
PSO parameters has been widely studied in the related 
literature (Clerc & Kennedy, 2002; Trelea, 2003). Based 
on this literature, the PSO parameters as shown in 
Table 8 were selected for this case study. Ten sets of 
optimization results were obtained by PSO. The 
standard deviation of the fitness of the ten 
independent runs is about 0.02. Table 9 shows the best 
fitness value of PSO out of the ten independent runs is 
- 33.7. 

The optimized set-points recommended by PSO are 
also denormalized into true values and reorganized 
into six compensation maps. The optimal maps are 
shown in Table 10. 

Evaluation of the Results 

By comparing the fitness value of PSO with the GA, 
the fitness value of PSO (-33.7) is larger than that of the 
GA (-39.5). Therefore, the overall compensation 
performance using the set-points produced by PSO is 
better than that based on the GA. By comparing the 
standard deviation of the ten independent runs of PSO 
with the GA, PSO is less than the GA. In other words, 
the PSO results are more stable than those with the GA 
because PSO is less sensitive to the initial values. 
Another issue is about the computational time. The 
computational times for MIO LS-SVCM + GA and 
MIO LS-SVCM + PSO are almost the same. On a 2.66 
GHz Intel Core 17-920 PC with 6 GB RAM on-board, 
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both methods spent approximately 76 hours for 
modelling and optimization. This result reveals that 
both computer-aided methods can greatly shorten the 
calibration time than the traditional manual approach 
that usually spends many months. Although the 
computational time of PSO is no advantage over the 
GA, PSO can provide the maximum fitness and the 
minimum standard deviation of the fitness. Hence, 
PSO is recommended to be the optimization tool for 
this problem domain. The aforementioned results 
show that PSO is superior to the GA in this application. 
The reasons for this superiority may be that the PSO 
mechanism generates a more diverse population 
during the whole iteration process. 

The actual output performance using the set-points 
recommended by PSO was also loaded to the test car 
and evaluation tests were performed using the chassis 
dyno. The overall accuracy of the optimal output 
performance 'Op' relative to actual test results 'Da' is 
94.2%, again verifying that the ECCS models using the 
MIO LS-SVCM method are accurate and reliable. The 
overall actual improvement of Da over Dbest, which is 
the result of the best fitness among the 300 sample 
datasets, is 12.96%. The results are promising and they 
are presented in Table 11. 

Conclusions 

ECMC is a very difficult problem so that very little 
research on computer-aided optimization for ECMC 
can be found in the existing literature. In this paper, a 
new framework called MIO LS-SVCM combining 
committee machine, MIO LS-SVM and evolutionary 
optimization for ECMC has been successfully 
developed and implemented. In the case of ECMC, the 
data consists of tunable and untunable parameters. It 
is unwise to include untunable parameters such as IAT, 
ET, BP, and BV for modelling and optimization. With 
reference to the idea of committee machine, this study 
firstly proposes to divide such complex system model 
into a number of simpler sub-models based on 
reasonable combinations of untunable parameters, so 
that the ECC sub-models do not contain any 
untunable parameter. With this division, the MIO LS- 
SVM in the framework can be employed to produce a 
set of ECC sub-models. Then, the models are 
systematically weighted and combined to form a 
multi-objective function for optimization. Finally, two 
common evolutionary optimization techniques, GA 
and PSO, are applied to the objective function to 
search for optimal values for calibration maps. 
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Evaluation tests show that the overall prediction 
accuracy of the LS-SVCM models is very good, at 
94.1%. Experimental results also show that PSO is 
superior to GA under the LS-SVCM. Based on the 
optimal calibration maps obtained from PSO, the 
overall system performance relative to Dbest is 
improved by 12.96%. Furthermore, the output 
performance using the recommended optimal maps is 
in good agreement with the engine real performance 
(about a 94.2% fit). Therefore, the proposed approach 
of MIO LS-SVCM + PSO is effective and performs well 
in this application. From the perspective of automotive 
engineering, this paper is the first in the literature that 
applied computational intelligence approach to solve 
the ECMC problem, and the proposed intelligent 
approach can greatly shorten the calibration time than 



the traditional manual approach. 

In addition, the research is the first attempt at 
integrating diverse paradigms (Latin-hypercube 
sampling, nonlinear regression, MIO LS-SVCM and 
GA/PSO) into a general framework for large-scale 
constrained multivariable optimization problems 
under the conditions of unknown system model and 
the existence of untunable and tunable system 
parameters. The successful integration of GA and PSO 
with MIO LS-SVCM also indicates that the proposed 
modelling and optimization framework can fit many 
evolutionary optimization techniques. It is believed 
that the proposed framework can be applied to other 
similar engineering modelling and optimization 
problems as well. 



TABLE 7 OPTIMAL COMPENSATION MAPS WITH THE GA 
Optimal Fuel Compensation Map for ET 



Eire 


10 


20 30 40 50 60 70 


80 


90 


100 
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33 


32 26 20 17 14 7 


3 





3 


5 


Optimal Ignition Compensation Map for ET 
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Optimal Fuel Compensation Map for IAT 
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Optimal Ignition Compensation Map for IAT 
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Optimal Fuel Compensation Map for BV 



BV/V 



11 



12 



13 



14 



Optimal Fuel Compensation Map for cranking 



Cranking time/ms 



0-200 



TABLE 8 PSO PARAMETERS 



Number of generation 


1000 


Population size 


50 


Inertial weight 


0.9 


Cognitive parameter 


2 


Social parameter 


2 
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TABLE 9 PSO RESULTS OF DIFFERENT OUTPUT FUNCTIONS 
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Engine output 


Knh(x)C/o) 


CTi(x)(s) 


AAi(x) 


AT<(x)(Nm) 


function Mh,i(x) 


h=l 


h=2 


h=3 


h=4 


Mh,20,20,n(x) 


15.77 


0.33 


0.08 


-0.08 


Mh,30,20,n(x) 


15.85 


3.06 


0.08 


-0.67 


Mh,40,20,12(x) 


16.04 


1.55 


0.07 


6.41 


MhM,20,n(x) 


15.89 


1.14 


0.05 


2.33 


Mh,60,20,12(x) 


15.90 


3.76 


0.05 


-1.66 


Mh,70,20,u(x) 


16.39 


0.51 


0.07 


0.68 


MhsomMx) 


16.00 


2.30 


0.05 


0.91 


Mh,60,30,12(x) 


16.00 


2.09 


0.06 


0.27 


Mh,70M,12(x) 


16.00 


2.81 


0.47 


0.60 


MhmmMx) 


15.77 


2.34 


0.48 


2.09 


Mh,90,30,12(x) 


16.35 


1.45 


0.79 


1.92 


Mh,10030,12(x) 


16.28 


0.80 


0.50 


0.27 


Mh,90,4B,12(x) 


16.39 


1.70 


0.06 


0.01 


Mh,W0A0,l2(x) 


15.85 


2.46 


0.96 


1.10 


Mh,100f0,12(x) 


16.22 


2.02 


0.15 


0.68 


Mh,20,20Mx) 


15.87 


1.13 


0.11 


0.11 


Mh,30,20,t3(x) 


15.95 


2.06 


0.11 


-0.47 


Mh,4B,20,13(x) 


16.14 


0.55 


0.01 


6.66 


Mhfi0,20,13(x) 


15.99 


1.14 


0.03 


2.53 


Mh,60,20,13(x) 


16.00 


1.76 


0.08 


-1.46 


Mh,70,20,13(x) 


16.49 


0.51 


0.10 


0.88 


Mh,50,30,13(x) 


16.10 


1.30 


0.07 


1.11 


Mh,60,30,13(x) 


16.10 


1.01 


0.09 


0.47 


Mh,70,30,13(x) 


16.10 


2.81 


0.50 


0.81 


Mh,S0,30,13(x) 


15.70 


0.79 


0.32 


1.89 


Mh,90,30,13(x) 


16.13 


1.56 


0.15 


2.12 


A/I u mn 3/1 7?fv 1 


16.38 


0.85 


0.53 


-1.46 


M 11,90,40,13 (x) 


16.50 


1.70 


0.09 


0.20 


Mh,W0A0,13(x) 


16.71 


1.19 


0.13 


2.27 


Mh,100M13(x) 


16.32 


2.02 


0.17 


0.88 


Overall average 


16.11 


1.62 


0.22 


1.05 


Fitness value 






-33.7 
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TABLE 10 OPTIMAL COMPENSATION MAPS WITH PSO 
Optimal Fuel Compensation Map for ET 



ET/'C 10 20 30 40 50 60 70 80 90 100 110 




Optimal Ignition Compensation Map for ET 

ET/*C 10 20 30 40 50 60 70 80 90 100 110 

(Ignition angle 3.5 3 2.5 2 1.5 1 0.5 -2 -4 

trim) 



Optimal Fuel Compensation Map for IAT 



IAT/°C 10 20 30 40 50 60 

(%fUelt " m) 7 5 2 -2 -5 " 

Optimal Ignition Compensation Map for IAT 



AT/*C 10 20 30 40 50 60 



Optimal Fuel Compensation Map for BV 



BV/V 11 


12 13 


14 


(|asec) 1124 1025 895 790 


Optimal Fuel Compensation Map for cranking 




Cranking time/ms 


0-200 










TABLE 11 COMPARISON BETWEEN OPTIMIZATION RESULTS, ACTUAL TEST RESULTS AND Dbest 


knk ( o /o) 


CT {S) AX 


a ryi Overall 
A I (Nm) 


Actual test results using the set-points ^ 
recommended by PSO, Da 


1.50 0.24 


1.12 


Optimization results using PSO, Op 16.11 


1.62 0.22 


1.05 


Optimization results using GA, O g 16.42 


1.75 0.52 


1.02 



Optimization results using the best 
set-points among the 300 sample 
datasets, Dbest 



18.41 



1.80 



0.28 



1.23 



Accuracy (Op relative to Da) 

99.38% 92.00% 91.67% 93.75% 94.20% 

|0-D a | 
(1 — - — - — ^)xl00% 

Actual improvement (Da relative to 

Dbest) 
D bes, ~ D a xlQ0 % 



11.95% 16.67% 14.29% 8.94% 12.96% 
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List of Abbreviation 

BP Barometric pressure 

BV Battery voltage 

CT Cranking time 

DOE Design of experiment 

ECC Engine compensation control 

ECCS Engine compensation control system 

ECMC Engine compensation map calibration 

ET Engine temperature 

EMS Engine management system 

GA Genetic algorithm 

IAT Inlet air temperature 

LS-SVM Least-squares support vector machine 

LS-SVCM Least-squares support vector committee 
machine 



MIO 

NN 

PSO 



Multi-input/output 

Neural network 

Particle swarm optimization 
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